Studies of spin-orbit scattering in noble-metal nanoparticles using energy level 
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The effects of spin-orbit scattering on discrete electronic energy levels are studied in copper, sil- 
ver, and gold nanoparticles. Level- to- level fluctuations of the effective ^-factor for Zeeman splitting 
are characterized, and the statistics are found to be well-described by random matrix theory predic- 
tions. The strength of spin-orbit scattering increases with atomic number and also varies between 
nanoparticles made of the same metal. The spin-orbit scattering rates in the nanoparticles are in 
order of magnitude agreement with bulk measurements on disordered samples. 

PACS numbers: 73.22.-f, 72.25.Rb, 73.23.Hk 
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Spin-orbit (SO) interactions are central to determin- 
ing the energy-level structure and dynamics of atomic 
nuclei, large-atomic-number atoms, and a wide variety 
of solid-state systems. They have important effects on 
electronic band structures jjj], they influence the sym- 
metry class of random-matrix theories applied to quan- 
tum systems Q , and they are put to use in proposals for 
controllable manipulation of electron spins || . A partic- 
ularly fundamental method for probing SO interactions 
in solid-state systems is to examine individual quantum 
energy levels in nanoparticles, both because SO matrix 
elements can be measured directly and because SO in- 
teractions produce strong modifications of the effective 
g-factors for spin-Zeeman splitting. Previously, measure- 
ments of three individual energy levels in particular gold 
nanoparticles showed that SO effects reduce the ^-factor 
to g ~ 0.28 — 0.45 in sharp contrast with Al par- 
ticles, which ordinarily give g ~ 2 j^]. Adding 4% gold 
to an aluminum nanoparticle also increased the SO scat- 
tering, thereby reducing the g-factors in the particle to 
^ ~ 0.5 - 0.8 §. In this paper we report systematic 
studies of SO effects in noble-metal nanoparticles cov- 
ering a range of atomic number (Z) - copper, silver, 
and gold. We characterize the mesoscopic fluctuations, 
present within single particles, of the g-factors and SO 
matrix elements (H QO ), and find that the statistics for 
the g-factors are well-described by recent predictions of 
random matrix theory (RMT) @,^|. From the values of 
the g-factors and the mean energy- level spacing 5 |Q, 
we extract SO scattering rates for each material, and 
we find order of magnitude agreement with expectations 
based on weak localization and electron spin resonance 
measurements. 

SO interactions cause the electronic energy levels in a 
metal grain to be not purely spin-up or spin-down, and as 
a consequence they reduce the g-factor for spin-Zeeman 
splitting below the free-electron value of 2. In low-Z 
nanoparticles for which SO scattering is sufficiently weak 
(with (H so ) <C S) perturbation theory can be used to 
calculate the g-factor for a given electronic state n [fol : 
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Here the \I/ m denote the unperturbed (pure-spin) single- 
electron eigenstates with energy E m . It is expected that 
the SO matrix elements will vary between different pairs 
of unperturbed energy levels because they depend on the 
precise nature of the wavefunctions, which are chaotic 
and fluctuating [|TT[] . The local energy level spacing will 
also vary. As a result, g n should exhibit mesoscopic fluc- 
tuations for different levels, n. For heavier elements, per- 
turbation theory ceases to be valid. However, for large 
SO-interaction strengths, Brouwer et al. and Matveev et 
al. have recently employed RMT to predict the statistical 
distribution P(g) of g-factors in a nanoparticle @,^]. In 
j?J, the strength of SO interactions is tuned continuously, 
allowing numerical studies of the full crossover from the 
Guassian orthogonal ensemble (GOE) to the Gaussian 
symplectic ensemble (GSE) for the energy-level statis- 
tics. In this crossover, it is found that the mean value of 
P(g) shifts to smaller g-factors and the variance of the 
distribution exhibits a maximum for (g) ~ 1.2. 

Our tunneling samples are fabricated in a vertical ge- 
ometry through a 3-10 nm bowl-shaped hole etched in 
a SiaN4 membrane filial (see inset, Fig. 1). One elec- 
trode is made by evaporating 1750 A of Al on top of 
the membrane to fill the bowl. Following oxidation in 
50 mTorr of O2 for 3 minutes to form a tunnel barrier, 
5-20 A of Cu, Ag, or Au is evaporated on the lower side 
of the device to form particles with radii of 3-5 nm. A 
second tunnel barrier is then formed by the evaporation 
of 11 A of AI2O3. This is followed by 1750 A of Al to 
make the lower tunneling electrode. The samples are 
screened at 4.2 K to select those whose current- volt age 
(I-V) curves exhibit Coulomb-staircase steps, indicating 
tunneling through a single nanoparticle. Then more de- 
tailed I-V curves are measured in a dilution refrigerator, 
with all electrical leads passing through copper-powder 
filters to minimize the effects of RF radiation. 

Figure 1 shows I-V and dl/dV-V curves for a Cu 
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FIG. 1. (a) I-V and (b) dl/dV-V curves for Cu #1. (c) 
Magnetic field dependence of peaks in dl/dV for energy levels 
exhibiting avoided level crossings. The V scale has been con- 
verted to energy, using the procedure discussed in the text. 
Inset: Cross-sectional device schematic. 

device (Cu#l). As V is increased, / increases in discrete 
steps due to tunneling through well-resolved energy lev- 
els within the nanoparticle ||. The resonances in the 
differential conductance directly correspond to the en- 
ergy level spectrum of the nanoparticle. In an applied 
magnetic field (H) most resonances in dl/dV divide into 
2 peaks, allowing measurements of the energy-difference 
for spin-Zeeman splitting between Kramers-pair eigen- 
states, AE = gjisH. In order to determine the effective 
g-factors, we must convert the magnitude of the splitting 
from the measured value in voltage to energy, correcting 
for capacitive division of the voltage across the two tun- 
nel junctions according to AE = AV/ (1 + C1/C2) where 
C1/C2 is the capacitance ratio for the two tunnel junc- 
tions. This ratio can be determined by measuring the 
shift in V for the resonances as the aluminum leads are 
changed from superconducting to normal in small mag- 
netic fields, or by comparing the voltage for tunneling 
through the same quantum state for positive and nega- 
tive V |J. The low- field g-factors for sample Cu#l fall 
in the range from 1.30±0.04 to 1.82±0.05. 

In the weak SO-scattering limit where perturbation 
theory is valid, the effects of SO interactions on the mag- 
netic field evolution of neighboring energy levels (shown 
in detail in Fig. 1(c) and in broader view in Fig. 2(a)) can 
be modeled by diagonalizing a simple 2x2 Hamiltonian 
matrix with off-diagonal matrix elements (H so ) coupling 
the eigenstates |6j. The minimum energy difference at 
the avoided level crossing directly gives 2(H so ) between 
those two states. The solid lines in Fig. 1(c) correspond 
to solutions of 2 x 2 Hamiltonians with (H so )=76 /ieV for 
the avoided crossing near 4.9 meV, and (iJ so )=134 /ieV 
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FIG. 2. Color scale plots of dl/dV in (a,b) Cu, (c) Ag, 
and (d) Au samples as a function of energy and magnetic field. 
The arrows denote the avoided crossings of Fig. 1. 

for the avoided crossing near 5.4 meV. These matrix el- 
ements are in the size range expected based on the mea- 
sured g-factors and the mean level spacing S = 0.70 meV. 
If the level spacing were uniform and equal to this value, 
and all the matrix elements (H so ) were ~ 130 /ieV, then 
Eq. (1) would predict g = 2 - 6.6((H so ) /5) 2 ~ 1.77. If 
energy-level fluctuations are included, the estimate for 
the typical g-factor is somewhat smaller ||. In accord 
with Eq. (1), we have seen that in a given sample the 
states which happen to have smaller- than- average energy 
spacings from neighboring levels generally have smaller- 
than-average g-factors. 

For the heavier elements Ag and Au, the g-factors be- 
come smaller with increasing atomic number. Figure 2 
shows color scale plots of dl/dV for increasing magnetic 
field in Cu, Ag, and Au samples. We find g-factors in the 
range 1.3 to 1.9 for Cu, 0.25 to 1.1 for Ag, and 0.05 to 0.19 
in Au for the samples shown. Table 1 lists the measured 
mean values of the g-factors and their standard devia- 
tions for all the samples we have measured with more 
than three resolvable resonances. 

The most direct way to compare the statistics of the 
measured g-factors with theoretical predictions is by 
plotting the integrated probability distribution of the g- 
factors for each sample. In Fig. 3 the points correspond 
to the g-factors from several Au, Ag, and Cu samples, 
while the solid lines are the theoretical predictions from 
Brouwer et al. with the SO-scattering strength pa- 
rameter A adjusted to minimize the least squares error 
between the theoretical and experimental distributions. 
Our results are in good qualitative agreement with the 
RMT predictions. In particular, for a single value of 
the adjustable parameter A, the theoretical distributions 
account reasonably for both the mean value of g and the 
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FIG. 3. Integrated probability distributions of the 
^-factors for several samples (points), compared to the pre- 
dictions of random-matrix theory with the SO-interaction 
strength adjusted for the best fit (lines). See Table 1 for 
the sample parameters. 

standard deviation for each sample. 

We do observe significant differences in the average 
value of the ^-factor for different nanoparticles composed 
of the same material, e.g., for the two Cu samples de- 
picted in Fig. 3 and for the Ag samples #1 and #2. This 
is true even though the mean level spacings for the two 
Cu samples are nearly identical. Because of these differ- 
ences, we cannot meaningfully consolidate the data from 
many samples to improve the quality of our statistics 
in making comparisons to the RMT predictions. The 
physical origin of SO interactions at low temperatures 
within noble metal nanoparticles is expected to be scat- 
tering from defects or the nanoparticle surface [|l3|-[l5|. 
Sample-to-sample variations in the strength of SO inter- 
actions are therefore reasonable in nanoparticles, if the 
defect concentration or the surface-to- volume ratio may 
vary, for instance due to size and shape differences be- 
tween particles. 

A few of the energy levels shown in Fig. 2 do not ex- 
hibit Zeeman splitting. We observe this only in samples 
where we can identify the tunneling transitions as cor- 
responding to changes from an odd to an even number 
of conduction electrons within the particle. We make 
this identification by observing the lowest energy tunnel- 
ing resonance. If the nanoparticle initially has an odd 
number of electrons, then the first energy level accessible 
for tunneling will have only one spin state available, and 
therefore the lowest energy transition does not undergo 
Zeeman splitting (Fig. 2(a,c,d)) ||. In other words, the 
tunneling transition occurs via a spin-singlet state. Sim- 
ilarly, we can explain the lack of observed splitting in 
some higher energy resonances as due to tunneling via 



excited spin-singlet states. (See the states near 7 meV in 
Fig. 2(c) and 0.3 meV in Fig. 2(d).) It is also possible 
that the Zeeman splitting may sometimes be sufficiently 
small for a level in Ag or Au that we cannot resolve it. 
Our resolution limit is approximately g > 0.03. From 
the theoretical probability distributions of Fig. 3, about 

0. 02% of the states in Ag and 1.4% of the states in Au 
can be expected to have g<0.03. 

In some samples, background-charge noise and 
nonequilibrium effects can interfere with accurate deter- 
minations of the statistics of the ^-factors. For instance, 
as a function of the voltage across the device or as a 
function of time, the background charge affecting the 
nanoparticle can sometimes take on an ensemble of dif- 
ferent values, shifting the energy levels that one is try- 
ing to measure. When this happens, current steps due 
to background-charge changes may be misinterpreted as 
due to quantum states, or the same quantum state may 
be measured several times at different values of voltage. 
In one Cu sample not shown in this paper, we measured 
six apparently different resonances all with an identical 
^-factor of 1.50, which we ascribe to this effect. This 
behavior was not seen in any other samples measured in 
this study. If several closely-spaced levels are observed 
in a sample, ways to diagnose the cause as background 
noise include the presence of negative differential conduc- 
tance even with normal-state electrodes, the presence of 
apparent level crossings rather than avoided crossings, 
and identical g-f actors for several resonance peaks. 

An important question for understanding the proper- 
ties of nm-scale devices is whether the strength of SO 
scattering is consistent with expectations based on mea- 
surements made by other techniques on larger samples, or 
whether there is new physics at small scales. Because ef- 
fective ^-factors measured by different experimental tech- 
niques can probe physically distinct quantities || , we will 
make the comparison based on estimates for SO scatter- 
ing times (t so ) that we can extract from our data. We 
will also make comparison only to samples with mean 
free paths for electron scattering comparable to the size 
of our nanoparticles, since this scattering is the probable 
source of SO-induced spin relaxation [[l3|Jl5|| . In the limit 
of strong SO scattering in a nanoparticle, the relationship 
between the ^-factor and r so has been solved analytically 

(g 2 ) = ( 2 ) 

For the Au samples, using the values of 5 listed in Table 

1, we find l/r so = (2-8) x lO 12 *" 1 . The Ag and Cu sam- 
ples fall outside the strong-scattering regime. However by 
fitting to the measured ^-factors the numerical distribu- 
tions calculated by Brouwer et al, we can still determine 
r so for each sample from the single fitting parameter A 
defined in j?]]: l/r so = \ 2 3/(ti?i). By this method, the 
SO scattering rates determined for the two Ag samples 
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are quite different, l/r so = 3 x 10 n s 1 and 2 x 10 12 s x , Nanofabrication Users Network, 
while for Cu l/r so = (2 - 8) x lO 11 ^ 1 . 

For Cu, we can also use perturbation theory to check 

the rate determined from the RMT approach. The value 
of t so within perturbation theory is Eg] 
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For (H so ) ~ 130/ieV and 6=0.7 meV, we find l/r so ~ 
2 x 10 n s _1 , in agreement with the rate for Cu#l listed 
in Table 1. 

These values for r so can be compared to the SO scatter- 
ing rates measured previously by weak localization, and 
from the linewidths of conduction electron spin resonance 
(CESR) in metal foils. From weak localization measure- 
ments on quench-condensed noble-metal films, Bergmann 
@ found 1/t S o = 5.4 x lO 12 *- 1 for Au, 2.6-3.1 x lO 11 *- 1 
for Ag (depending on the disorder) and 3.2 x 10 n s -1 for 
Cu. The values for Au and Cu agree to within a factor 
of three with all the rates we extract, while the weak- 
localization rate for Ag coincides with sample Ag#2 but 
not Ag#l. The CESR comparisons are less direct, be- 
cause these measurements are performed on crystalline 
foils with thicknesses on the order of microns or more. 
In these experiments it is generally found that the spin- 
flip rate scales inversely with sample thickness due to the 
dominant effect of surface scattering |0| . The measured 
values include spin-flip rates of 7 x Krs -1 for Au foils 
one hundred microns thick ||l9| , 1.5 x 10 8 s _1 for a 5-/zm- 
thick Ag foil H, and 2.6 x lO 8 ^ 1 for a 1.2 /im Cu film 
|^of . Assuming that the rates would scale as 1 / thickness 
down to a thickness of 5 nm, we arrive at estimates of 
1 x lO 14 ^- 1 for Au, 1 x lO 11 ^- 1 for Ag, and 6 x lO 10 *" 1 for 
Cu. These SO scattering rates are also in rough order-of- 
magnitude agreement with the values obtained from the 
^-factor analysis. 

In conclusion, we have measured the ^-factor distribu- 
tions in Cu, Ag, and Au nanoparticles using tunneling 
spectroscopy. The ^-factors within a given sample ex- 
hibit mesoscopic fluctuations, with distributions in good 
accord with random-matrix-theory predictions. We have 
observed significant differences in both the ^-factor distri- 
butions and r so between different nanoparticles composed 
of the same metal, illustrating the importance of the spe- 
cific sample structure in the origin of SO scattering. The 
SO scattering rates we extract for the nanometer-scale 
particles are in rough order-of-magnitude agreement with 
extrapolations from previous CESR and weak localiza- 
tion measurements. 
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TABLE I. Sample parameters. N is the number of reso- 
nances resolved, (g) is the mean ^-factor, cr e (g) is the experi- 
mental standard deviation, a t (g) is the standard deviation of 
the theory curve with the best-fit value of the SO parameter 
A as denned in [7], 6 is the mean level spacing, and l/r so is 
the spin-orbit scattering rate calculated as discussed in the 
text. 
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